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TRANSVERSE  KELVIN-HELMHOLTZ  INSTABILITY  WITH  PARALLEL 
ELECTRON  DYNAMICS  AND  COULOMB  COLLISIONS 


I.  INTRODUCTION 

Velocity  shear  layers  with  steep  gradient  scale  lengths  are 
frequently  observed  in  the  high  latitude  auroral  zone  (Kelley  and 
Carlson,  1982;  3asu  et  al.,  1985;  Kelley  and  Earle,  1986).  Recent  HILAT 
satellite  measurements  reported  by  Basu  et  al.  (1985)  show  flows  with 
moderate  to  strong  velocity  shears  with  a  shallow  irregularity  spectrum; 
often  field  aligned  currents  are  also  observed  in  the  regions  of  sheared 
velocity  flows.  The  observations  suggest  that  the  east-west  E  x  B  drift 


velocity, 

for 

example,  is 

inhomogeneous  and 

usually  reverses 

i  uO 

direction 

as 

one  moves  in 

the 

north-south  dir 

ection.  This  kind 

of 

velocity 

shear 

transverse  to 

the 

magnetic  field 

is  usually  a  source 

of 

Kelvin-Helmholtz  instability  (Chandrasekhar ,  1961;  Mikhailovski i , 
1972).  In  the  conventional  Kelvin-Helmholtz  instability  the  mixing  of 
fluids  with  different  velocities  is  confined  to  two  dimensions 
transverse  to  the  magnetic  field;  the  ion  polarization  drift  leads  to 
the  instability  while  the  electrons  simply  E  *  3  drift  in  the 
perpendicular  plane.  In  the  auroral  ionosphere  one  cannot  simply  ignore 
the  dynamics  parallel  to  the  magnetic  field  because  there  could  be 
additional  sources  of  free  energy  or  damping  that  may  affect  the 
transverse  mixing  of  the  fluids  with  sheared  velocity.  Our  aim  is  to 
introduce  the  third  dimension  into  the  two  dimensional  analysis  of 
Kelvin-Helmholtz  instability  and,  as  a  first  step,  we  study  the  effects 
of  parallel  electron  dynamics.  Furthermore,  since  the  electron  parallel 
motion  couples  the  collisionless  and  collisions!  regimes  of  the 
ionosphere,  we  also  include  Coulomb  collisions  between  the  electrons  and 
ions.  A  preliminary  analysis,  in  the  collisionless  domain,  was  performed 
by  Thompson  (1983)  who  showed  that  the  parallel  electron  motion  could 
have  stabilizing  influence  on  the  Kelvin-Helmholtz  instability.  In  this 
paper  we  perform  f  i  >  a  more  complete  analysis  and  present  a  general 
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stability  criterion  for  a  specific  velocity  profile,  (ii)  present  the 
effects  of  collisions  on  the  instability.  In  analogy  with  the  role  of 
gravity  on  Kelvin-Helmholtz  instability  in  neutral  fluids,  we  observe 
that  the  parallel  electron  dynamics  simulate  buoyancy  of  a  fluid  in  a 
gravitational  field. 

II.  THEORY 

We  consider  a  homogeneous  plasma  that  is  drifting  across  a  magnetic 

A 

field  (B  ■  B  z)  because  of  an  equilibrium  inhomogeneous  electric  field 

A 

(E  -  Eq(x)  x).  We  choose  EQ(x)  =  Eg  tanh(x/L)  where  L  is  the  scale 
length  of  the  shear  layer.  In  this  paper  we  restrict  the  analysis  to 
the  domain  L  >>  p^  and  j<w|  <<  ,  where  and  n  are  the  ion 
gyroradius  and  gyrofrequency ,  respectively.  Incompressible  plasma 
motion  in  two  dimensions  across  an  uniform  magnetic  field  is  determined 
by  the  constraint  that  the  perpendicular  current  be  divergence  free, 


:i ) 


where  ^  »  ne(v^  -  v  ) ,  e  is  the  charge,  n  is  the  plasma  density,  v 

and  v^  are  the  perpendicular  drifts  of  electrons  and  ions.  The  veloci¬ 
ties  vg^  and  are  obtained  from  the  momentum  equations  by  assuming 
the  electrons  only  E  *  B  drift,  while  the  ions  experience  both  the 
E  *  B  drift  and  the  polarization  drift.  Linearizing  (1)  with  a 

A  A 

perturbed  electrostatic  potential  of  the  form  $  -  <p ( x )  exp  [-i(wt  “ 

k  y ) ]  with  frequency  u>  and  wavenumber  k  yields  the  equation  governing 
y  / 

the  Kelvin-Helmholtz  i.nstabilty: 
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where  Vq(x)  »  -  (c/B)Eg(x).  This  is  a  well  studied  equation 

(Chandrasekhar ,  1961;  Mikhailovskii ,  1961*;  Cuzdar  et  al . ,  1982; 
Satyanarayana  et  al . ,  1983);  instability  (Y  >  0,  where  Y  is  given 
by  uj  «  +  iY)  occurs  in  the  wavenumber  domain  0  <  k^L  <  1  for  Vg(x)  = 

V q  tanh(x/L) . 

To  extend  the  above  theory  to  three  dimensions,  we  need  to  consider 


V- 

CN 


m 

■w 


which  allows  for  a  parallel  current  along  the  field  lines.  In  this 
paper  we  consider  perturbations  along  the  magnetic  field  and  also 
equilibrium  currents  along  the  field  lines.  We  include  electron 
collisions  and  use  a  BGX  model  to  calculate  the  electron  susceptibility 
(Clemmow  and  Dougherty,  1969).  Kinetic  treatment  of  the  electron 
parallel  motion  yields  a  general  equation  for  the  perturbed  electro- 

A  A 

static  potential  $  -  q(x)  exp[-i(wt  -  k  z  -  k  y ) ]  in  the  domain 

2  y 

j ui |  and  Vq/L  «  n. 
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where  =  Qga^,  is  t!le  electron  (ion)  gyro-frequency ,  ve  i3  the 

electron  thermal  velocity,  £  ■  (u  -  k  V  -  k  V  ♦  iv  )/k  v  ,  ,  is 

e  y  0  zd  e;  z  e  d 

the  equilibrium  electron  drift  along  the  magnetic  field,  and  va  is  the 
electron  collision  frequency.  For  va  *  0,  (M)  is  the  same  as  equation 
(8)  of  Thompson  (1983)  without  the  ion  terms;  we  can  ignore  the  parallel 
ion  terms  if  we  restrict  our  analysis  to  p^/L  <  1  where  is  the  ion 
gyroradius.  We  have  introduced  the  parallel  electron  drift,  Vd,  and 
collisional  effects  which  are  absent  in  Thompson  (1983). 


III.  COLLISIONLESS  REGIME 

(a)  Analytical  Results 

In  this  section  we  set  v 


0.  A  simple  stability  boundary 


defining  the  regions  of  k  and  k  where  the  system  is  stable  or  unstable 

Z  j 

can  be  easily  obtained  in  the  domain  where  the  electrons  behave  as  a 
fluid;  this  is  the  region  where  the  argument  of  the  Z  function  is  large 


(i.e. ,  IS, 


>>  1 


-  1/5.  -  1/2  5e 
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In  this  domain  the  Z  function  can  be  approximated  as 
With  this  approximation  (U)  takes  the  form 


.  I  ilia— 

3x2  y  (uJ  ky/0)  2  (u>  -  kyVQ) 


?]♦  -  0. 


o  / 


Equation  (5)  is  identical  to  the  Rayleigh  equation  governing  the 

interchange  of  a  drifting,  weakly  inhomogeneous  fluid  in  a  gravitational 

field  with  a  sheared  drift  velocity  (Drazin,  1953;  Zhandrasexhar , 

2  o  o 

1961).  Here  we  identify  the  term  w  <”/ k^  with  the  buoyancy  term  that 

'6n  t*  y 

is  driving  the  Rayleign-Taylor  instability  g/Lp,  where  g  is  gravity  and 
Ln  is  the  density  gradient  scale  length.  In  the  spirit  of  the  neutral 
fluid  proolem,  we  define  a  dimensionless  number 


J 


(6) 
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2 


,2 
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2 
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Once  we  have  made  this  identification,  it  follows  that  for  VQ(x) 
tanh(x/L)  the  stability  boundary  is  simply  (Drazin,  1958) 


J 


k2L2(l 

y 


(7) 


The  modes  are  stable  for  J  >  1/4  or 


K  ,  Vn  2 

—  >  -  f— —) 

k2  2  W  ‘ 

y 


(8) 


Thus,  parallel  electron  dynamics  has  a  stabilizing  influence  on  the 
Kelvin-Helmholtz  instability. 


Figure  1  shows  the  stability  boundary  where  we  plot  J  versus  k  L. 
We  note  that  the  conclusions  of  Thompson  (1983)  are  qualitatively 
correct.  However,  our  analysis  shows  that  J  <  1/4  for  instability  as 
opposed  to  J  <  1  given  by  Thompson  (1983)  who  did  not  consider  the 
explicit  shear  profile  V(x)  »  VQ  tanh(x/L).  A  derivation  of  this 
stability  criterion  based  upon  the  energy  principle  is  given  in  Sec. 
III.c. 


(b)  Numerical  Results 

We  now  solve  (4)  numerically  to  further  show  the  effects  of 

the  parallel  dynamics.  In  Figure  2  we  plot  the  normalized  growth 

2  2  2 

rate  y/(vql]  vs.  a  (*  a»^L  /v*)  for  k^L  »  0.5  (which  gives  maximum 

growth  for  k2  »  0),  kz/ky  *  1.0  x  1 0~ ,  Vg/ve  =  0.01,  and  «  0.  Curve 
A  shows  the  growth  rate  for  the  collisionless  case,  while  Curve  S  shows 
the  growth  rate  when  electron  collisions  are  included  (to  be  discussed 
in  Sec.  IV).  For  a  »  0  we  find  Y/(v  /L )  =  0.19  which  is  the  usual 
result.  As  at  is  increased,  we  see  the  Kelvin-Helmholtz  growth  rate 

♦I 

drops.  Since  for  large  ot  the  Kelvin-Helmholtz  term  (i.e.  «  '/„)  does  not 

o 

contribute  to  (4),  the  mode  equation  essentially  becomes 


;  v, 


which  allows  only  stable  roots  in  the  fluid  limit  with  u  -  k..V~  = 

r  s  u 

(1/2)(k  /k  )uu  .  From  Curve  A  of  Figure  2  we  also  see  that  the  modes 
z  y  in 

are  stable  for  a  >  9.  This  is  in  agreement  with  the  criterion  J  >  1  /4 
previously  discussed.  We  find  that  for  reasonable  values  of  the 
parallel  drift  velocity  (i.e.,  Vd  <  v  )  these  results  are  not 

significantly  altered.  We  find  that  in  this  regime  the  argument  of  the 

Z  function  is  large,  indicating  that  the  electrons  are  merely  acting  as 

a  fluid  and  wave  particle  effects  are  not  playing  an  important  role.  On 

the  other  hand,  for  sufficiently  small  a  the  theory  breaks  down  as 

shown  by  the  kink  in  Fig.  2;  the  orbits  begin  to  play  an  important  role 
here  (L  -  p^)  and  kinetic  effects  need  to  included. 

We  now  solve  (4)  for  several  mode  numbers.  Figure  3  shows  the 

normalized  growth  rate  Y/(VgL)  as  a  function  of  k  L  for  VQ/ve  =  0.01, 

a  ■  1 ,  and  k^L  =  0,  1.0  x  10'3,  2.0  x  10‘3,  3.0  x  10"3,  and  4.0  x  10“3 

as  curves  A,  B,  C,  D  and  E,  respectively.  Curve  A  shows  the 

conventional  Kelvin-Helmholtz  mode;  Y  >  0  in  the  wavenumber  domain  0  < 

k,,L  <  1  and  Y  »  0.19  for  k„L  *  0.46.  Several  points  are  to  be  noted 
y  max  y 

from  this  figure:  (i)  as  k  L  is  increased  the  wavenumber  of  the 

Z  _o 

maximally  growing  mode  increases;  for  example,  for  kzL  =  4.0  x  10  the 

growth  rate  maximizes  at  kyL  ~  0.8,  (ii)  the  growth  steadily  decreases 

as  k  L  increases,  dropping  by  50$  for  k  L  =  3  x  10~3,  and  (iii)  there  is 
"  z 

a  narrower  range  of  k  L  for  unstable  modes;  for  k_L  =  4.0  x  10~3  the 

y  £ 

waves  are  unstable  only  in  the  region  0.60  <  kyL  <  0.93. 

Equation  (4)  is  solved  both  for  the  eigenvalues  and  the 
eigenfunctions  (i.e.,  the  perturbed  electrostatic  potentials).  The 
eigenfunctions  reveal  some  features  of  the  stabilizing  influence  of  the 
electron  parallel  motion.  Briefly,  the  method  used  to  solve  (4)  is  the 
following;  The  asymptotic  form  of  the  solution  is  assumed  to  be  WK3 
type 

~  TTu  exp  [“  !  ^  dxl 

Q 

where  Q  is  the  coefficient  of  4  in  (4).  For  a  given  value  of  k  L,  we 

a  j 

search  for  an  eigenvalue  ui  such  that  integration  of  (4)  with  41  as  the 

X 

boundary  condition  yields  a  localized  solution  whose  logarithmic 


derivative  3/3x(!n  $)  is  continuous  across,  for  example,  the  origin.  We 


note  that  the  electron  terms  are  competing  with  -  k^  in  (4)  at  x  =  »; 
if  the  electron  terms  are  dominating  Q  at  x  =  «  then  Q  could  become 
positive,  leading  to  oscillatory  type  solution  at  x  =  «.  On  the  other 


hand,  if  the  electron  terms  are  not  large  such  that  Q  at 


=  »  is 


negative,  that  is 


ky,  the  WKB  solutions  are  exponentially  decaying. 


Figure  4  shows  the  solutions  of  (4)  for  k  L  =  3*0  x  10~^  and  k„L  = 

z  / 

0.7  for  which  Y/fv^/L)  =  0.09  and  w  /fv./L.)  =  0.  We  see  that  both  the 

(j  r*  u 

real  (A)  and  imaginary  (B)  parts  are  exponentially  decaying  at  x 

*  +  »  and  are  fairly  localized  near  the  origin  where  the  velocity 

reversal  is  taking  place.  The  localization  width  is  proportional  to  the 

thickness  of  the  shear  layer,  Ax  -  4L.  The  value  of  Q  at  x  =  ®  is 

2  2 

negative  but  greater  than  -  k  L  .  In  Figure  4b  we  show  the  full 

y  •» 

perturbed  electrostatic  potential,  $  -  $(x)  exp  [ik  y]  using  the  wave 

*  J 

function,  <f>(x),  for  the  maximally  growing  mode  of  Figure  4a.  The  dashed 
(solid)  lines  indicate  the  relative  negative  (positive)  values  of  $.  We 
note  that  these  values  correspond  to  positive  and  negative  vorticities 
of  the  fluid;  the  phase  of  $  is  such  that  it  allows  mixing  of  the  fluid 

with  different  vorticities  thereby  causing  the  perturbation  to  grow. 


In  Figure  5  we  show  the  eigenfunctions  for  the  lower  cutoff 

wavenumber  k  L  a  0.51  corresponding  to  curve  D  of  Fig.  3  for  which  k„L  = 

J  *-• 

3.0  x  10-3;  the  growth  rate  and  the  real  frequency  for  this  case  are 
Y/(Vq/l)  =  0  and  ^/(V^/L)  =  -  0.057,  respectively.  In  this  case  we 

find  that  Q  at  x  =  -«  is  positive  showing  oscillatory  type  solutions, 
while  Q  is  negative  at  x  =  *  giving  rise  to  an  exponential Ly  decaying 
solution.  In  comparison  with  Figure  4a,  this  figure  shows  that 
electrons  are  convecting  energy  away  from  the  localized  region  in  the 
shear  layer,  thus,  leading  to  stablization  of  the  instability;  this 

point  is  discussed  in  detail  in  the  next  section.  In  Figure  5b  we  show 

the  corresponding  $.  Comparison  of  Figs.  4b  and  5b  shows  that  the 
parallel  motion  of  electrons  changes  the  phase  of  t>  and  the  finite  real 
frequency  introduces  a  relative  velocity  to  the  wave  preventing  any 
mixing  of  the  fluid  in  the  x  <  0  and  x  >  0  regions,  thus  causing  the 
perturbation  not  to  grow. 


.  mVV-V. 
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c.  Energy  Principle 

In  this  section  we  present  energy  principle  arguments  to  show 
the  effects  of  electron  parallel  dynamics.  The  Kel vin-Helmholtz 
instability  arises  basically  due  to  the  conversion  of  available  kinetic 
energy  of  relative  motion  of  the  equilibrium  flow  into  wave  energy.  In 
this  case  the  electrons  are  confined  to  the  plane  perpendicular  to  the 
magnetic  field  and  only  E  *  3  drift.  If  we  allow  the  electrons  to  move 
along  the  magnetic  field,  the  available  kinetic  energy  in  the 

equilibrium  flow  has  to  overcome  the  parallel  electron  compressional 
energy  for  instability  to  occur. 

We  can  quantify  these  arguments  and  obtain  a  stability  condition  as 
follows.  We  consider  the  interchange  of  two  neighboring  volume  elements 
at  x  and  x  +  5x  which  are  moving  with  velocities  V0  and  V q 

+  5V,  respectively.  The  available  kinetic  energy  is 

5T  =»  ~  nm.[v2  ♦  (VQ  +  5V)2  -  \  (2VQ  *  SV)2]  =  \  nm.[5V)2  (10) 

where  m^  is  the  mass  of  the  ions  and  n  is  the  density.  The  compres¬ 
sional  energy  associated  with  fluid  moving  a  distance  Sz  (i.e.,  parallel 

to  B)  is  given  by 

<5z 

5W  =  <  I  nm  5v  dz  >  =  ~  nm  <5v  5z>  (11) 

‘  e  z  2  e  z 

where  m0  is  the  electron  mass  and 

6v  -  -  (e/m  )3p/3z  (12) 

z  e 

and  <  >  denotes  time  averaged  quantities.  Since  the  perturbed 

electron  motion  in  the  perpendicular  plane  is  predominantly  E  *  3 
motion,  we  have 


5v  =  -  (c/B)3$/3y  . 


From  (12)  and  (13)  wo  have 


and 


dv 


•  k  Q 
dv  z  e 
z 


_  r> 

dx  _  y 

dv  k  n 
z  z  e 


Substituting  (14)  in  (11)  we  obtain 


i  k_  8.  i  ><8 

<5W  =  -  nma  <dz  dvx>  =  ^  nm^  (-2-2)  <dx  dvz> 

. y  y 


from  which  the  compressional  energy  is  given  as 

i  k  8^  -5  o 

5W  =  —  nm  (-f^)2  (dx)2  . 


(14) 


(15) 


(16) 


i  nmjdv)2  <  1  nme(^i)2  (fix}1 


From  (10)  and  (16)  the  stability  condition,  dT  <  6W,  is 

k  a 

Z  £ 
c 

y 

or 

(?r)2  <  2  (Vi 


(17) 


(18) 


since  dV 


(aVQ/3x)  dx. 


Thus  , 


the  stability  condition  can  be  written  as 


> 


1_ 

2 


(19) 


which  agrees  with  (8)  derived  in  Section  II. 
IV.  C0LLISI0NAL  REGIME 


In  this  section  we  study  the  effects  of  electron  collisions  or.  the 
transverse  Kelvin-Heimhoitz  instability.  In  Fig.  6  we  plot  the 
normalized  growth  rate  (y/(V ,/U)  as  a  function  of  k„L  for  y„  =  , 

v0/ve  3  G.01,  k,L  =  1.0  x  10  3,  and  k^L  =  0.6.  The  variation  of  toe 
normalized  growth  rate  for  y  =  0  is  shown  by  curve  9 ,  and  for 
'■'o  *  100  (V^/L)  by  curve  3.  We  note  that  the  electron  collisions 
introduce  a  rew  scale  size,  the  collisior.al  mean  free  path 


8 


.-..a 


-■.v, 


I 


A  _  (=  v  /v  ).  Since  A  ,  <<  A  ,  electron  motion  parallel  to  3  is 

mfp  e  e  tnfp  z 


inhibited,  and  as  a  consequence  the  k^L  domain  over  which  the  modes  are 
how  unstable  is  much  larger.  This  is  shown  clearly  in  Fig.  6,  where  the 
collisional  growth  rate  (curve  B)  goes  to  zero  for  kzL  >  0.07,  while  the 
collisionless  growth  rate  (curve  A)  drops  to  zero  sharply  for  kzL  > 
0.005. 


In  addition,  curve  B  of  Fig.  2  shows  the  collosional  growth  rate 
versus  the  coupling  parameter  a.  For  the  range  of  a  considered, 
electron  collisions  make  the  growth  rate  almost  independent  of  a.  The 
electron  collisions  thus  have  a  destabilizing  influence  on  Kelvin- 
Helmholtz  instability.  This  would  in  turn  alter  the  stability  boundary 
in  such  a  way  that  the  basic  Kelvin-Helmholtz  mode  is  unstable  in  a 
larger  k_  domain  than  the  collisionless  case. 


V.  SUMMARY  AND  DISCUSSION 

In  this  paper  we  have  considered  the  effect  of  parallel  electron 
motion  on  the  transverse  Kelvin-Helmholtz  instability  as  a  first  step 


toward  a  three-dimensional  analysis  of  the  nature  of  convective  flows  in 


the  high  latitude  ionosphere.  In  the  collisionless  regime  (va  =  0),  we 
show  that  for  an  inhomogeneous  velocity  profile  given  by  V(x)  =  tanh 


(x/L)  the  Kelvin-Helmholtz  instability  is  stable  for  Kz/k-/ 


>  ( V  /Lai. J  k  l(  1  -  k2L2)1/2  in  the  limit  a  -  k  V.  >>  k7va.  We  further 
o  Jtn  y  y  y  u  c.  - 


show  that  the  physical  mechanism  for  stabilization  is  the  compressional 


energy  given  to  the  electrons  parallel  to  3.  In  the  collisional  regime 


(  v„  *  0),  we  show  that  the  Kelvin-Helmholtz  instability  is  not  as 

t  ^ 


easily  stabilized  by  finite  k  effects;  this  is  attributed 


inhibition  of  parallel  electron  motion  f or  A  „  <<  A 

mf  p  ; 


We  apply  this  stability  criterion  for  parameters  reieva 
high  latitude  ionosphere.  We  take  w 


ass um t 


k  /'<  =*  (o  /<j  J 


1  /  2 


Ih 


3.2  x  1 0d  sec"1 ,  kyL  =  0. 


arc 


where  and  are 
parallel  conductivities,  respectively  (Farley,  1959).  Th 


2  y  l  \ 


he  perpend! cular  and 
va  i  ue 


of  a  /a.  varies  with  altitude  in  the  ionosphere  and  is  generally  taker, 

-L  I  _C 

to  be  in  the  range  10  ■'  (lower  F  region)  -  id  '  '  ;  re 


region , 
- 1  . 


10  '  v  t  o  ps 

Thus,  for  instability  to  occur  it  is  required  that  V^/L  >2.3  sec 
the  topside  region.  Velocity  shears  have  seen  reported  in  the  rang 
Vq/L  =  0.5  -  20  Hz  (Kelley  and  Carlson,  1977;  Earle  and  Kelley,  1 986'  s 


tnat 


ir.sfaoii  itv 


DOUi  ' 


.  y  b- 


act  i  ve 


;os  i.  >- 


Ir.  the  iow^r  F  region,  ^here  v  -  10"  -  ’ 


and  V „  '  v . 

j  - 


10  sec  are  needed  to  drive  the 


10 


-2 


velocity  shears  of  order  1  - 


Kelvin-Helmholtz  modes  unstable  (e.g.,  for  k^L  =  0.6  and  kzL  =  0.01). 

However,  these  conclusions  are  predicated  on  the  important  assumption 
.  ,  i  /o 

that  k,/kv  »  o  /o,  ~  which  may  not  be  the  case.  Recent  analyses  of 

z  y  1  | 

unstable  flute  modes  in  barium  clouds  (i.e.,  the  gradient  drift 


instability)  (Sperling  et  al.,  1984;  Drake  et  al . ,  1  985)  indicate  that 

k_/k„  =*  0  within  the  unstable  cloud.  A  more  detailed  three-dimensional 
z  y 

analysis  (Drake  and  Huba,  1987)  is  needed  to  assess  the  role  of  the 
Kelvin-Helmholtz  instability  in  high  latitude  dynamics.  Finally,  we  add 
that  we  have  neglected  ion-neutral  collisions  in  the  analysis;  these  can 
be  important  in  the  lower  F  region.  In  a  forthcoming  3tudy  (Mitchell  et 
al . ,  1987)  we  will  present  analytical  and  computational  results 

demonstrating  the  effect  of  ion-neutral  collisions  on  the  transverse 
Kelvin-Helmholtz  instability.  In  the  linear  regime,  it  is  found  that 
ion-neutral  collisions  have  a  stabilizing  influence  on  the  instability. 
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stability  boundary  in  Che  fluid  limit, 
marginal  stability  boundary  (y  =  0)  of  J 
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imaginary  parts  of  <f>,  respectively,  (b)  Contour  plot 


Fig.  5)  Eigenfunction  for  the  mode  corresponding  to  Curve  C  of 
for  k  L  »  O.uS  (y  *  0).  (a)  Curves  A  and  3  refer  to  ti 

J  A 

and  imaginary  parts  of  j,  respectively.  Contour 

of  $  -  $(x)  exp  [ik  y 3. 
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ligcnfunction  for  the  mode  corresponding  to  Curve  I)  of  Pig. 
A,  /.  “  0  48  (y  —  0).  (a)  Curves  A  and  B  refer  to  the  real 
imaginary  parts  of  <£,  respectively,  (b)  Contour  plots 
>  exp  |/A,.iy|. 
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01CY  ATTN  CODE  4  6 1 
0  1  C  Y  ATTN  CODE  402 
01CY  ATTN  CODE  420 
01CY  ATTN  CODE  421 

COMMANDER 

AEROSPACE  DEFENSE  C  0 MM AN  D / X P  D 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFS,  CO  80912 
01CY  ATTN  XPDQQ 
0 1 C Y  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 


AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-P ATTERSON  AFB ,  OH  45433 
01CY  ATTN  AAD  WADE  HUNT 
01 CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  &  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
0 1 CY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 

u  a  u  cr  nu  *c*d  u  •  *▼->«  « <r  ss  o 

01CY  ATTN  J.  DEAS 
ESD/SCD-4 

COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  45433 
01CY  ATTN  N I  CD  LIBRARY 

01CY  ATTN  ETDP  B.  BALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB  .  NY  13  4m 

01CY  ATTN  DOC  LIBRARY/TSLD 
01CY  ATTN  OCSE  V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AFB ,  NB  68  1  1  3 
0 1 CY  ATTN  XPFS 


HANSCOM 

AFB, 

MA  01731 

SAMSO/MN 

01  CY 

ATTN 

OPR 

HAROLD 

GARDNER 

NORTON  AF  3 ,  CA  92409 

01  CY 

ATTN 

LK8 

(  MINUTEMAN  ) 

KENNETH  S 

,W.  CHAMPION 

0 1 C  Y  ATTN  MNNL 

01  CY 

ATTN 

OPR 

ALVA  T 

.  STAIR 

01  CY 

ATTN 

PHD 

JURGEN 

BUCHAU 

COMMANDER 

01  CY 

ATTN 

PHD 

JOHN  P 

.  MULLEN 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 

HANSCOM  AFB,  MA  01731 

AF  WEAPONS  LABORA 

TORY 

0 1 CY  ATTN  EE  P  A.  LORENTZEN 

KIRTLAND 

AFT, 

NM 

87117 

01  CY 

ATTN 

SUL 

DEPARTMENT  OF  ENERGY 

01  CY 

ATTN 

CA 

ARTHUR 

H.  GUENTHER 

LIBRARY  ROOM  G-042 

01  CY 

ATTN 

NT  YC 

E  1  LT  . 

G.  KRAJEI 

WASHINGTON,  DC  20545 

0 1 C  Y  ATTN  DOC  CON  FOR  A.  LABOWI 

AFTAC 

PATRICK 

AFB  , 

FL  3 

2925 

DEPARTMENT  OF  ENERGY 

0  1  CY 

ATTN 

TN 

ALBUQUERQUE  OPERATIONS  OFFICE 

P  .  0  .  SOX  5  4  00 
ALBUQUERQUE ,  NM  37115 

01CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 
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EG&G ,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  855KH 

0  1  C  Y  ATTN  DOC  CON  FOR  J.  3RESDL0VE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  9«550 

01  CY  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

0 1 C Y  ATTN  DOC  CON  FOR  L-389  R.  OTT 

0 1 C Y  ATTN  DOC  CON  FOR  L-j!  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  875^5 

01CY  ATTN  DOC  CON  FOR  J.  WOLCOTT 

01 CY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 

0 1 C  Y  ATTN  DOC  CON  FOR  E.  JONES 

01CY  ATTN  DOC  CON  FOR  J.  MALIK 

0 1 C  Y  ATTN  DOC  CON  FOR  R.  JEFFRIES 

01 CY  ATTN  DOC  CON  FOR  J.  2INN 

01CY  ATTN  DOC  CON  FOR  D.  WESTERVELT 

01CY  ATTN  D.  SAPP  ENFIELD 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D«33 

LOS  ALAMOS.  NM  3 7 5 « 5 
01CY  ATTN  S.P.  GARY 
01CY  ATTN  J.  BOROVSKY 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  205^5 

01  CY  ATTN  DOC  CON  DR.  YO  SONG 

NATIONAL  OCEANIC  A  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

0 1 C  Y  ATTN  R.  GRUB3 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 


P.O.  BOX 

9295 

7 

LOS  ANGELES, 

CA 

90009 

01  CY 

ATTN 

T  _ 

GARFUNKEL 

01  CY 

ATTN 

7  _ 

SALMI 

01  CY 

ATTN 

V  . 

JOSEPHSON 

0  1  CY 

ATTN 

s . 

BOWER 

01  CY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CO  HP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN.  TX  78758 

01CY  ATTN  L.  SLOAN 
01 CY  ATTN  R.  THOMPSON. 


SANDIA  LABORATORIES 
P.O.  BOX  5800 


.BUQUERQUE, 

NM 

37115 

01  CY 

ATTN 

DOC 

CON 

FOR 

W.  3R0WN 

01  CY 

ATTN  DOC  CON 
THORNBROUGH 

FOR 

A  . 

01  CY 

ATTN 

DOC 

CON 

FOR 

T.  WRIGHT 

01  CY 

ATTN 

DOC 

CON 

FOR 

D.  DAHLGREN 

01  CY 

ATTN 

DOC 

CON 

FOR 

3  l  a  1 

01  CY 

ATTN 

D  I  V 

DOC 

CON 

FOR 

SPACE  PROJE 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  }U550 

01  CY  ATTN  DOC  CON  FOR  B.  MURPHEY 
0 1 C  Y  ATTN  OOC  CON  FOR  T.  COOK 


BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  3  OX  983 
BERKELEY.  CA  9  U  7  0 1 


01  CY 

ATTN 

J  . 

,  WORKMAN 

01  CY 

ATTN 

f> 

PRETTIE 

01  CY 

ATTN 

r* 

»J  . 

.  BRECHT 

BOEING 

COMPAN 

V 

THE 

P.O.  BOX  3707 
SEATTLE,  WA  98 1 

2P 

0  1  CY 

ATTN 

c. 

KEISTER 

01  CY 

ATTN 

D  . 

MURRAY 

0  1  CY 

ATTN 

G. 

HALL 

01  CY 

ATTN 

J  . 

KENNEY 

CHARLES  STARK  DRAPER  LABORATORY.  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
01 CY  ATTN  D. B  COX 
0 1 C  Y  ATTN  J.P.  GILMORE 


COMSAT  LABOR ATOR IES 
22300  COMSAT  DRIVE 
CLARKSBURG.  MD  20871 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING. 
ITHACA,  NY  14850 

0 1 C  Y  ATTN  D.T.  FARLEY.  JR. 

ELECTROSPACE  SYSTEMS.  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
0 1 C Y  ATTN  H.  LOGSTON 
0 1 C  Y  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wllanire  3 1  vd . 

Santa  Monica,  CA  90401 
01CY  ATTN  C.B.  GABBARD 
0 1 C  Y  ATTN  R.  LELEVIER 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  3555 
PHILADELPHIA,  PA  19101 
0 1 C  Y  ATTN  M.H.  30RTNER 
SPACE  SCI  LAS 

GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

0 1 C  Y  ATTN  G.  MILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

0 1 C Y  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 C  Y  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

0 1 C Y  ATTN  DICK  STEINHOF 


HSS ,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  6 1 3  2  0 

(ALL  CORRES  ATTN  DAN  MCCLELLAND 
0 1  C  Y  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  V  A  22311 
0 1 C  Y  ATTN  J.M.  AEIN 
0 1 C Y  ATTN  ERNEST  3AUER 
0 1 C  Y  ATTN  HANS  WOLFARD 
0 1 C Y  ATTN  JOEL  BENGSTON 

I NTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY  ,  NJ  07110 

01 CY  ATTN  TECHNICAL  LIBRARY 

J AYCOR 

11011  TORREY ANA  ROAD 
P.O.  BOX  351  5^ 

SAN  DIEGO,  CA  921 38 

0 1 C  Y  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  _IBRARIAN 
01CY  ATTN  THOMAS  POTE.MRA 
01 CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEACHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDI ES 

316  STATE  STREET  (P.O  DRAWER  QQ ) 
SANTA  3ARSARA,  CA  93102 
01CY  ATTN  D  AS  I  AC 
0 1 C  Y  ATTN  WARREN  3.  KNAPP 
0 1 C  Y  ATTN  WILLIAM  MCNAMARA 
0 1 C  Y  ATTN  B.  GAMBILL 


LINKABIT  CORP 
t  0 « 5 3  ROSELLE 
SAN  DIEGO ,  CA  92121 

0  1  C  Y  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC 
P.O.  BOX  50« 

SUNNYVALE,  CA  9^083 
0  1  C Y  ATTN  DEPT  60-12 
0 1 C  Y  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  4  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  9*4  3014 

0 1 C Y  ATTN  MARTIN  WALT  DEPT  52-12 
0 1 C Y  ATTN  W.L.  IMHOF  DEPT  52-12 
0 1  C  Y  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

0 1  C Y  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32305 

0 1 C Y  ATTN  R.  HEFFNER 

MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  3EACH,  CA  926*47 
0 1 C  Y  ATTN  N.  HARRIS 
01 CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MR02 
01CY  ATTN  W  .  OLSON 
0 1  C  Y  ATTN  R.W.  HALPRIN 
0 1 C Y  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
0 1 C  Y  ATTN  P.  FISCHER 
0 1 C Y  ATTN  W.F.  CREVIER 
01 CY  ATTN  STEVEN  L.  CUTSCHE 
0 1 C  Y  ATTN  R.  30GUSCH 
0 1 C  Y  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KIL3 
0 1 C Y  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
0 1 C  Y  ATTN  CONRAD  L .  LONGMIRE 
0 1 C  Y  ATTN  5.  WHITE 
01CY  ATTN  R.  STAGAT 


MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM  87106 
0 1 C Y  R.  STELLINGWERF 
0 1 C Y  M.  ALME 
0  1  C Y  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01 CY  ATTN  W.  HALL 
0 1 C Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1  2  3  **  0  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

0 1 C  Y  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

0 1 C  Y  ATTN  IONOSPHERIC  RESEARCH  LA 

PHOTOMETRICS,  INC. 

9  ARROW  drive 
WOBURN,  MA  01801 

0 1 C Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01 CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  9 U 6 1 0 
ATTN  A.  THOMSON 

R  4  D  ASSOCIATES 

P.O.  30 X  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 

01CY  ATTN  WILLIAM  B.  WRIGHT,  JR. 

01CY  ATTN  WILLIAM  J.  KAR2AS 

01CY  ATTN  H.  ORY 

01CY  ATTN  C.  MACDONALD 

01CY  ATTN  BRIAN  LAMB 

0 1 C  Y  ATTN  MORGAN  GROVER 


RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY ,  MA  01 776 

0 1 C  Y  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK.  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  INCORPORATED 
1 1 50  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

01CY  ATTN  LEWIS  M.  LINSON 
0 1 C  Y  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
0 1 C Y  ATTN  CURTIS  A.  SMITH 


TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 


REDONDO 

BEACH , 

CA  90278 

OlCY 

ATTN 

R  . 

K.  PLEBUCH 

OlCY 

ATTN 

S. 

ALTSCHULER 

01  CY 

ATTN 

D. 

DEE 

OlCY 

ATTN 

D/ 

STOCKWELL 

SNTF / 1 575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01303 
0 1 C  Y  ATTN  W.  REIDY 
OlCY  ATTN  J.  CARPENTER 
0 1 C  Y  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

OlCY  ATTN:  N.  ZABUSKY 


SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
OlCY  J.  COCKAYNE 
OlCY  E.  HYMAN 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 


MENLO 

PARK,  ' 

CA 

94025 

OlCY 

ATTN 

J  . 

CASPER 

OlCY 

ATTN 

DONALD  NEILSON 

OlCY 

ATTN 

ALAN  BURNS 

OlCY 

ATTN 

G. 

SMITH 

OlCY 

ATTN 

R. 

TSUNODA 
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ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  C.  CHESNUT 

01  CY 

ATTN 

CHARLES  L.  RINO 

OlCY 

ATTN 

WALTER  JAYE 

OlCY 

ATTN 

J  . 

VICKREY 

OlCY 

ATTN 

RAY  L.  LEADABRAND 

01  CY 

ATTN 

G  . 

CARPENTER 

01  CY 

ATTN 

G. 

PRICE 

OlCY 

ATTN 

R  . 

LIVINGSTON 

01  CY 

ATTN 

V  . 

GONZALES 

OlCY 

ATTN 

D. 

MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OlCY  ATTN  W.P.  80QUIST 
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